The role of body size and cuticular hydrocarbons in the desiccation resistance of invasive Argentine ants (Linepithema humile).

An insect's cuticle is typically covered in a layer of wax prominently featuring various hydrocarbons involved in desiccation resistance and chemical communication. In Argentine ants (Linepithema humile), cuticular hydrocarbons (CHCs) communicate colony identity, but also provide waterproofing necessary to survive dry conditions. Theory suggests different CHC compound classes have functional trade-offs, such that selection for compounds used in communication would compromise waterproofing, and vice versa. We sampled sites of invasive L. humile populations from across California to test whether CHC differences between them can explain differences in their desiccation survival. We hypothesized that CHCs whose abundance was correlated with environmental factors would determine survival during desiccation, but our regression analysis did not support this hypothesis. Interestingly, we found the abundance of most CHCs had a negative correlation with survival, regardless of compound class. We suggest that the CHC differences between L. humile nests in California are insufficient to explain their differential survival against desiccation, and that body mass is a better predictor of desiccation survival at this scale of comparison.

Reactive oxygen species-mediated immunity against bacterial infection in the gut of cat fleas (Ctenocephalides felis).

Fleas (Order Siphonaptera) transmit numerous bacterial pathogens that cause severe human diseases (e.g., cat scratch disease, flea-borne spotted fever, murine typhus, plague). Because initial entry of these infectious agents occurs while blood feeding, the immune response in the flea gut is considered to be the first line of defense against invading microbes. However, relatively few studies have identified the flea immune molecules that effectively resist or limit infection in the gut. In other hematophagous insects, an immediate immune response to imbibed pathogens is the generation of reactive oxygen species (ROS). In this study, we utilized cat fleas (Ctenocephalides felis) to investigate whether oral infection with a well-known insect bacterial pathogen (Serratia marcescens) induces ROS synthesis in the flea gut, and whether production of ROS provides a defense mechanism against microbial colonization. Specifically, we treated fleas with an antioxidant to limit the number of free radicals in the digestive tract prior to infection, and then measured the following: S. marcescens infection loads, hydrogen peroxide (ROS) levels, and mRNA abundance of ROS signaling pathway genes. Overall, our data shows that ROS levels increase in response to infection in the flea gut, and that this increase helps to strengthen the flea immune response through the microbicidal activity of ROS.

Lack of parent-of-origin effects in Nasonia jewel wasp: A replication and extension study.

In diploid cells, the paternal and maternal alleles are, on average, equally expressed. There are exceptions from this: a small number of genes express the maternal or paternal allele copy exclusively. This phenomenon, known as genomic imprinting, is common among eutherian mammals and some plant species; however, genomic imprinting in species with haplodiploid sex determination is not well characterized. Previous work reported no parent-of-origin effects in the hybrids of closely related haplodiploid Nasonia vitripennis and Nasonia giraulti jewel wasps, suggesting a lack of epigenetic reprogramming during embryogenesis in these species. Here, we replicate the gene expression dataset and observations using different individuals and sequencing technology, as well as reproduce these findings using the previously published RNA sequence data following our data analysis strategy. The major difference from the previous dataset is that they used an introgression strain as one of the parents and we found several loci that resisted introgression in that strain. Our results from both datasets demonstrate a species-of-origin effect, rather than a parent-of-origin effect. We present a reproducible workflow that others may use for replicating the results. Overall, we reproduced the original report of no parent-of-origin effects in the haplodiploid Nasonia using the original data with our new processing and analysis pipeline and replicated these results with our newly generated data.

Functional insights from the GC-poor genomes of two aphid parasitoids, Aphidius ervi and Lysiphlebus fabarum.

BACKGROUND: Parasitoid wasps have fascinating life cycles and play an important role in trophic networks, yet little is known about their genome content and function. Parasitoids that infect aphids are an important group with the potential for biological control. Their success depends on adapting to develop inside aphids and overcoming both host aphid defenses and their protective endosymbionts. RESULTS: We present the de novo genome assemblies, detailed annotation, and comparative analysis of two closely related parasitoid wasps that target pest aphids: Aphidius ervi and Lysiphlebus fabarum (Hymenoptera: Braconidae: Aphidiinae). The genomes are small (139 and 141 Mbp) and the most AT-rich reported thus far for any arthropod (GC content: 25.8 and 23.8%). This nucleotide bias is accompanied by skewed codon usage and is stronger in genes with adult-biased expression. AT-richness may be the consequence of reduced genome size, a near absence of DNA methylation, and energy efficiency. We identify missing desaturase genes, whose absence may underlie mimicry in the cuticular hydrocarbon profile of L. fabarum. We highlight key gene groups including those underlying venom composition, chemosensory perception, and sex determination, as well as potential losses in immune pathway genes. CONCLUSIONS: These findings are of fundamental interest for insect evolution and biological control applications. They provide a strong foundation for further functional studies into coevolution between parasitoids and their hosts. Both genomes are available at https://bipaa.genouest.org.

Sawfly Genomes Reveal Evolutionary Acquisitions That Fostered the Mega-Radiation of Parasitoid and Eusocial Hymenoptera.

The tremendous diversity of Hymenoptera is commonly attributed to the evolution of parasitoidism in the last common ancestor of parasitoid sawflies (Orussidae) and wasp-waisted Hymenoptera (Apocrita). However, Apocrita and Orussidae differ dramatically in their species richness, indicating that the diversification of Apocrita was promoted by additional traits. These traits have remained elusive due to a paucity of sawfly genome sequences, in particular those of parasitoid sawflies. Here, we present comparative analyses of draft genomes of the primarily phytophagous sawfly Athalia rosae and the parasitoid sawfly Orussus abietinus. Our analyses revealed that the ancestral hymenopteran genome exhibited traits that were previously considered unique to eusocial Apocrita (e.g., low transposable element content and activity) and a wider gene repertoire than previously thought (e.g., genes for CO2 detection). Moreover, we discovered that Apocrita evolved a significantly larger array of odorant receptors than sawflies, which could be relevant to the remarkable diversification of Apocrita by enabling efficient detection and reliable identification of hosts.

Desiccation Resistance and Micro-Climate Adaptation: Cuticular Hydrocarbon Signatures of Different Argentine Ant Supercolonies Across California.

Cuticular hydrocarbons (CHCs), the dominant fraction of the insects' epicuticle and the primary barrier to desiccation, form the basis for a wide range of chemical signaling systems. In eusocial insects, CHCs are key mediators of nestmate recognition, and colony identity appears to be maintained through a uniform CHC profile. In the unicolonial Argentine ant Linepithema humile, an unparalleled invasive expansion has led to vast supercolonies whose nestmates can still recognize each other across thousands of miles. CHC profiles are expected to display considerable variation as they adapt to fundamentally differing environmental conditions across the Argentine ant's expanded range, yet this variation would largely conflict with the vastly extended nestmate recognition based on CHC uniformity. To shed light on these seemingly contradictory selective pressures, we attempt to decipher which CHC classes enable adaptation to such a wide array of environmental conditions and contrast them with the overall CHC profile uniformity postulated to maintain nestmate recognition. n-Alkanes and n-alkenes showed the largest adaptability to environmental conditions most closely associated with desiccation, pointing at their function for water-proofing. Trimethyl alkanes, on the other hand, were reduced in environments associated with higher desiccation stress. However, CHC patterns correlated with environmental conditions were largely overriden when taking overall CHC variation across the expanded range of L. humile into account, resulting in conserved colony-specific CHC signatures. This delivers intriguing insights into the hierarchy of CHC functionality integrating both adaptation to a wide array of different climatic conditions and the maintenance of a universally accepted chemical profile.

The complete mitochondrial genome of the invasive Africanized Honey Bee, Apis mellifera scutellata (Insecta: Hymenoptera: Apidae).

The complete mitochondrial genome from an Africanized honey bee population (AHB, derived from Apis mellifera scutellata) was assembled and analyzed. The mitogenome is 16,411 bp long and contains the same gene repertoire and gene order as the European honey bee (13 protein coding genes, 22 tRNA genes and 2 rRNA genes). ND4 appears to use an alternate start codon and the long rRNA gene is 48 bp shorter in AHB due to a deletion in a terminal AT dinucleotide repeat. The dihydrouracil arm is missing from tRNA-Ser (AGN) and tRNA-Glu is missing the TV loop. The A + T content is comparable to the European honey bee (84.7%), which increases to 95% for the 3rd position in the protein coding genes.

Biased Allele Expression and Aggression in Hybrid Honeybees may be Influenced by Inappropriate Nuclear-Cytoplasmic Signaling.

Hybrid effects are often exhibited asymmetrically between reciprocal families. One way this could happen is if silencing of one parent's allele occurs in one lineage but not the other, which could affect the phenotypes of the hybrids asymmetrically by silencing that allele in only one of the hybrid families. We have previously tested for allele-specific expression biases in hybrids of European and Africanized honeybees and we found that there was an asymmetric overabundance of genes showing a maternal bias in the family with a European mother. Here, we further analyze allelic bias in these hybrids to ascertain whether they may underlie previously described asymmetries in metabolism and aggression in similar hybrid families and we speculate on what mechanisms may produce this biased allele usage. We find that there are over 500 genes that have some form of biased allele usage and over 200 of these are biased toward the maternal allele but only in the family with European maternity, mirroring the pattern observed for aggression and metabolic rate. This asymmetrically biased set is enriched for genes in loci associated with aggressive behavior and also for mitochondrial-localizing proteins. It contains many genes that play important roles in metabolic regulation. Moreover we find genes relating to the piwi-interacting RNA (piRNA) pathway, which is involved in chromatin modifications and epigenetic regulation and may help explain the mechanism underlying this asymmetric allele use. Based on these findings and previous work investigating aggression and metabolism in bees, we propose a novel hypothesis; that the asymmetric pattern of biased allele usage in these hybrids is a result of inappropriate use of piRNA-mediated nuclear-cytoplasmic signaling that is normally used to modulate aggression in honeybees. This is the first report of widespread asymmetric effects on allelic expression in hybrids and may represent a novel mechanism for gene regulation.

A Search for Parent-of-Origin Effects on Honey Bee Gene Expression.

Parent-specific gene expression (PSGE) is little known outside of mammals and plants. PSGE occurs when the expression level of a gene depends on whether an allele was inherited from the mother or the father. Kin selection theory predicts that there should be extensive PSGE in social insects because social insect parents can gain inclusive fitness benefits by silencing parental alleles in female offspring. We searched for evidence of PSGE in honey bees using transcriptomes from reciprocal crosses between European and Africanized strains. We found 46 transcripts with significant parent-of-origin effects on gene expression, many of which overexpressed the maternal allele. Interestingly, we also found a large proportion of genes showing a bias toward maternal alleles in only one of the reciprocal crosses. These results indicate that PSGE may occur in social insects. The nonreciprocal effects could be largely driven by hybrid incompatibility between these strains. Future work will help to determine if these are indeed parent-of-origin effects that can modulate inclusive fitness benefits.

Behavioural and genetic analyses of Nasonia shed light on the evolution of sex pheromones.

Sex pheromones play a pivotal role in the communication of many sexually reproducing organisms. Accordingly, speciation is often accompanied by pheromone diversification enabling proper mate finding and recognition. Current theory implies that chemical signals are under stabilizing selection by the receivers who thereby maintain the integrity of the signals. How the tremendous diversity of sex pheromones seen today evolved is poorly understood. Here we unravel the genetics of a newly evolved pheromone phenotype in wasps and present results from behavioural experiments indicating how the evolution of a new pheromone component occurred in an established sender-receiver system. We show that male Nasonia vitripennis evolved an additional pheromone compound differing only in its stereochemistry from a pre-existing one. Comparative behavioural studies show that conspecific females responded neutrally to the new pheromone phenotype when it evolved. Genetic mapping and gene knockdown show that a cluster of three closely linked genes accounts for the ability to produce this new pheromone phenotype. Our data suggest that new pheromone compounds can persist in a sender's population, without being selected against by the receiver and without the receiver having a pre-existing preference for the new pheromone phenotype, by initially remaining unperceived. Our results thus contribute valuable new insights into the evolutionary mechanisms underlying the diversification of sex pheromones. Furthermore, they indicate that the genetic basis of new pheromone compounds can be simple, allowing them to persist long enough in a population for receivers to evolve chemosensory adaptations for their exploitation.

Draft genome of the globally widespread and invasive Argentine ant (Linepithema humile).

Ants are some of the most abundant and familiar animals on Earth, and they play vital roles in most terrestrial ecosystems. Although all ants are eusocial, and display a variety of complex and fascinating behaviors, few genomic resources exist for them. Here, we report the draft genome sequence of a particularly widespread and well-studied species, the invasive Argentine ant (Linepithema humile), which was accomplished using a combination of 454 (Roche) and Illumina sequencing and community-based funding rather than federal grant support. Manual annotation of >1,000 genes from a variety of different gene families and functional classes reveals unique features of the Argentine ant's biology, as well as similarities to Apis mellifera and Nasonia vitripennis. Distinctive features of the Argentine ant genome include remarkable expansions of gustatory (116 genes) and odorant receptors (367 genes), an abundance of cytochrome P450 genes (>110), lineage-specific expansions of yellow/major royal jelly proteins and desaturases, and complete CpG DNA methylation and RNAi toolkits. The Argentine ant genome contains fewer immune genes than Drosophila and Tribolium, which may reflect the prominent role played by behavioral and chemical suppression of pathogens. Analysis of the ratio of observed to expected CpG nucleotides for genes in the reproductive development and apoptosis pathways suggests higher levels of methylation than in the genome overall. The resources provided by this genome sequence will offer an abundance of tools for researchers seeking to illuminate the fascinating biology of this emerging model organism.

The genome sequence of the leaf-cutter ant Atta cephalotes reveals insights into its obligate symbiotic lifestyle.

Leaf-cutter ants are one of the most important herbivorous insects in the Neotropics, harvesting vast quantities of fresh leaf material. The ants use leaves to cultivate a fungus that serves as the colony's primary food source. This obligate ant-fungus mutualism is one of the few occurrences of farming by non-humans and likely facilitated the formation of their massive colonies. Mature leaf-cutter ant colonies contain millions of workers ranging in size from small garden tenders to large soldiers, resulting in one of the most complex polymorphic caste systems within ants. To begin uncovering the genomic underpinnings of this system, we sequenced the genome of Atta cephalotes using 454 pyrosequencing. One prediction from this ant's lifestyle is that it has undergone genetic modifications that reflect its obligate dependence on the fungus for nutrients. Analysis of this genome sequence is consistent with this hypothesis, as we find evidence for reductions in genes related to nutrient acquisition. These include extensive reductions in serine proteases (which are likely unnecessary because proteolysis is not a primary mechanism used to process nutrients obtained from the fungus), a loss of genes involved in arginine biosynthesis (suggesting that this amino acid is obtained from the fungus), and the absence of a hexamerin (which sequesters amino acids during larval development in other insects). Following recent reports of genome sequences from other insects that engage in symbioses with beneficial microbes, the A. cephalotes genome provides new insights into the symbiotic lifestyle of this ant and advances our understanding of host-microbe symbioses.

Draft genome of the red harvester ant Pogonomyrmex barbatus.

We report the draft genome sequence of the red harvester ant, Pogonomyrmex barbatus. The genome was sequenced using 454 pyrosequencing, and the current assembly and annotation were completed in less than 1 y. Analyses of conserved gene groups (more than 1,200 manually annotated genes to date) suggest a high-quality assembly and annotation comparable to recently sequenced insect genomes using Sanger sequencing. The red harvester ant is a model for studying reproductive division of labor, phenotypic plasticity, and sociogenomics. Although the genome of P. barbatus is similar to other sequenced hymenopterans (Apis mellifera and Nasonia vitripennis) in GC content and compositional organization, and possesses a complete CpG methylation toolkit, its predicted genomic CpG content differs markedly from the other hymenopterans. Gene networks involved in generating key differences between the queen and worker castes (e.g., wings and ovaries) show signatures of increased methylation and suggest that ants and bees may have independently co-opted the same gene regulatory mechanisms for reproductive division of labor. Gene family expansions (e.g., 344 functional odorant receptors) and pseudogene accumulation in chemoreception and P450 genes compared with A. mellifera and N. vitripennis are consistent with major life-history changes during the adaptive radiation of Pogonomyrmex spp., perhaps in parallel with the development of the North American deserts.

Functional characterization of transcription factor motifs using cross-species comparison across large evolutionary distances.

We address the problem of finding statistically significant associations between cis-regulatory motifs and functional gene sets, in order to understand the biological roles of transcription factors. We develop a computational framework for this task, whose features include a new statistical score for motif scanning, the use of different scores for predicting targets of different motifs, and new ways to deal with redundancies among significant motif-function associations. This framework is applied to the recently sequenced genome of the jewel wasp, Nasonia vitripennis, making use of the existing knowledge of motifs and gene annotations in another insect genome, that of the fruitfly. The framework uses cross-species comparison to improve the specificity of its predictions, and does so without relying upon non-coding sequence alignment. It is therefore well suited for comparative genomics across large evolutionary divergences, where existing alignment-based methods are not applicable. We also apply the framework to find motifs associated with socially regulated gene sets in the honeybee, Apis mellifera, using comparisons with Nasonia, a solitary species, to identify honeybee-specific associations.

Recombination and its impact on the genome of the haplodiploid parasitoid wasp Nasonia.

Homologous meiotic recombination occurs in most sexually reproducing organisms, yet its evolutionary advantages are elusive. Previous research explored recombination in the honeybee, a eusocial hymenopteran with an exceptionally high genome-wide recombination rate. A comparable study in a non-social member of the Hymenoptera that would disentangle the impact of sociality from Hymenoptera-specific features such as haplodiploidy on the evolution of the high genome-wide recombination rate in social Hymenoptera is missing. Utilizing single-nucleotide polymorphisms (SNPs) between two Nasonia parasitoid wasp genomes, we developed a SNP genotyping microarray to infer a high-density linkage map for Nasonia. The map comprises 1,255 markers with an average distance of 0.3 cM. The mapped markers enabled us to arrange 265 scaffolds of the Nasonia genome assembly 1.0 on the linkage map, representing 63.6% of the assembled N. vitripennis genome. We estimated a genome-wide recombination rate of 1.4-1.5 cM/Mb for Nasonia, which is less than one tenth of the rate reported for the honeybee. The local recombination rate in Nasonia is positively correlated with the distance to the center of the linkage groups, GC content, and the proportion of simple repeats. In contrast to the honeybee genome, gene density in the parasitoid wasp genome is positively associated with the recombination rate; regions of low recombination are characterized by fewer genes with larger introns and by a greater distance between genes. Finally, we found that genes in regions of the genome with a low recombination frequency tend to have a higher ratio of non-synonymous to synonymous substitutions, likely due to the accumulation of slightly deleterious non-synonymous substitutions. These findings are consistent with the hypothesis that recombination reduces interference between linked sites and thereby facilitates adaptive evolution and the purging of deleterious mutations. Our results imply that the genomes of haplodiploid and of diploid higher eukaryotes do not differ systematically in their recombination rates and associated parameters.

Dll3 and Notch1 genetic interactions model axial segmental and craniofacial malformations of human birth defects.

Mutations in the Notch1 receptor and delta-like 3 (Dll3) ligand cause global disruptions in axial segmental patterning. Genetic interactions between members of the notch pathway have previously been shown to cause patterning defects not observed in single gene disruptions. We examined Dll3-Notch1 compound mouse mutants to screen for potential gene interactions. While mice heterozygous at either locus appeared normal, 30% of Dll3-Notch1 double heterozygous animals exhibited localized, segmental anomalies similar to human congenital vertebral defects. Unexpectedly, double heterozygous mice also displayed statistically significant reduction of mandibular height and decreased length of the [corrected] maxillary hard palate. Examination of somite-stage embryos and perinatal anatomy and histology did not reveal any organ defects, so we used microarray-based analysis of Dll3 and Notch1 mutant embryos to identify gene targets that may be involved in notch-regulated segmental or craniofacial development. Thus, Dll3-Notch1 double heterozygous mice model human congenital scoliosis and craniofacial disorders.

Identification of oscillatory genes in somitogenesis from functional genomic analysis of a human mesenchymal stem cell model.

During somitogenesis, oscillatory expression of genes in the notch and wnt signaling pathways plays a key role in regulating segmentation. These oscillations in expression levels are elements of a species-specific developmental mechanism. To date, the periodicity and components of the human clock remain unstudied. Here we show that a human mesenchymal stem/stromal cell (MSC) model can be induced to display oscillatory gene expression. We observed that the known cycling gene HES1 oscillated with a 5 h period consistent with available data on the rate of somitogenesis in humans. We also observed cycling of Hes1 expression in mouse C2C12 myoblasts with a period of 2 h, consistent with previous in vitro and embryonic studies. Furthermore, we used microarray and quantitative PCR (Q-PCR) analysis to identify additional genes that display oscillatory expression both in vitro and in mouse embryos. We confirmed oscillatory expression of the notch pathway gene Maml3 and the wnt pathway gene Nkd2 by whole mount in situ hybridization analysis and Q-PCR. Expression patterns of these genes were disrupted in Wnt3a(tm1Amc) mutants but not in Dll3(pu) mutants. Our results demonstrate that human and mouse in vitro models can recapitulate oscillatory expression observed in embryo and that a number of genes in multiple developmental pathways display dynamic expression in vitro.